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We simulate cluster formation of model colloidal particles interacting via DLVO (Derjaguin, Lan- 
dau, Vervey, Overbeek) potentials. The interaction potentials can be related to experimental con- 
ditions, denned by the pH-value, the salt concentration and the volume fraction of solid particles 
suspended in water. The system shows different structural properties for different conditions, includ- 
ing cluster formation, a glass-like repulsive structure, or a liquid suspension. Since many simulations 
are needed to explore the whole parameter space, when investigating the properties of the suspen- 
sion depending on the experimental conditions, we have developed a steering approach to control a 
running simulation and to detect interesting transitions from one region in the configuration space 
to another. The advantages of the steering approach and the restrictions of its applicability due to 
physical constraints are illustrated by several example cases. 



I. INTRODUCTION 

Soft matter physics is a large field which has 
gained more and more in importance during the 
last years. It comprises for example complex 
fluids, biological systems like membranes, so- 
lutions of large molecules like proteins, or sus- 
pensions of small soft or solid particles, which 
are commonly called "colloids" . Since one can 
find examples for these materials nearly every- 
where in everyday life (in medicine, food indus- 
try, paintings, glue, blood, ceramics,. . . ), many 
results in research have a high relevance for ap- 
plications. In this context, especially colloid 
science is one of the most active research fields 
of the present time. A considerable effort has 
been invested to describe colloidal suspensions 
from a theoretical point of view and by simula- 
tions [U |2l HJ [5] as well as to understand the 
particle-particle interactions [HI IZl EJ El EH EQ 
and the phase behavior [H HI El ED3 US]- 
Typical for soft matter physics, and especially 
for colloids, is that the important effects which 
govern the material properties are on a meso- 
scopic length scale, i.e., on lengths much larger 
than the atomistic scale, and much smaller than 
the macroscopic length scale. 

This mesoscopic level brings in the possibil- 
ity to control the parameters of a certain mate- 
rial by manipulating the processes on the meso- 
scopic length scale. For colloids there are many 
ways to control the interactions between the in- 
dividual colloidal particles. As an example, by 
changing the particle interactions, it is possi- 
ble to initiate a clustering process in a colloidal 
suspension [P71 \TE \ |H) 1 I5TJ 1 121 ]. |2"2 1 |2l] . 

However, on the mesoscopic length scale, of- 
ten many different effects are in a subtle in- 
terplay which makes it difficult to give quan- 
titative predictions. Computer simulations can 
help to study these systems in detail, e.g., to see 
the response of the system on a change in the 



particles' interactions. Such a change can be 
achieved by adding salt to a suspension or by 
changing the pH- value [HI QIJ [32j [24]. How- 
ever, in mesoscopic systems where many dif- 
ferent effects contribute to the overall behavior 
the parameter space is quite large and one has 
to perform a large number of different simula- 
tions, each of them requiring a lot of computing 
resources, to gain an understanding of all inter- 
relationships. 

Therefore, it is useful to be able to steer 
the simulation on-line. Steering in this context 
can mean to change the interaction potentials 
or other parameters like an externally applied 
shear rate, but it can also mean to go back in 
the simulation time and follow a different path 
starting from an earlier configuration, both in- 
duced by interaction with the user [231 125] . 

For our study we have modeled an aqueous 
suspension of Al 2 C"3-particles. The particles 
in our simulations are monodisperse spheres of 
0.37/im in diameter. They interact via DLVO 
(Derjaguin, Landau, Vervey, Overbeek) poten- 
tials [51 [7] as well as a repulsive force ensuring 
excluded volume for the particles. The prop- 
erties of this model system have been investi- 
gated by simulations and experiments in pre- 
vious works [H [Ml [27] OH HO W\- In this 
paper we discuss how steering of a simulation 
can help to explore the parameter space and 
which problems may occur when using steering 
techniques. 

In the following section we describe our sim- 
ulation method and focus on the implementa- 
tion of the steering. Then, we discuss some ex- 
amples of simulations where steering helps to 
explore the parameter space, but we also give 
examples in which a steered simulation might 
lead to different results if compared to an un- 
steercd simulation. Finally, we summarize our 
results and draw a conclusion. 
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II. SIMULATION METHOD 

Our simulation method is described in detail 
in Refs. [231 H3 121] and consists of two parts: a 
Molecular Dynamics (MD) code, which treats 
the colloidal particles, and a Stochastic Rota- 
tion Dynamics (SRD) simulation for the fluid 
solvent. In the MD part we include effective 
electrostatic interactions and van der Waals at- 
traction, known as DLVO potentials [6, 7 . The 
repulsive term results from the surface charge 
of the suspended particles 
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where d denotes the particle diameter, r the 
distance between the particle centers, e the el- 
ementary charge, T the temperature, &b the 
Boltzmann constant, and z is the valency of 
the ions of added salt. Sq is the permittivity of 
the vacuum, e r = 81 the relative dielectric con- 
stant of the solvent, n the inverse Debye length 
defined by k 2 = 8tt£bI, with ionic strength / 
and Bjerrum length £b — 7 A. We have related 
the effective surface potential C to the pH- value 
and the ionic strength of the solvent by means 
of a charge regulation model in our previous 
work [24 . Thus, the particle interaction po- 
tentials can be related to distinct experimental 
conditions. The second therm of the DLVO po- 
tentials which does not depend on the pH- value 
or the ionic strength is the attractive van der 
Waals interaction (A H = 4.76 • 10" 20 J is the 
Hamaker constant) [17] 
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The attractive contribution VVdw competes 
with the repulsive term and is responsible for 
the cluster formation one can observe for con- 
ditions in which the attraction dominates. 

Since DLVO theory is based on the assump- 
tion of large particle separations, it does not 
correctly reproduce the primary minimum in 
the potential, which should appear at particle 
contact. Therefore, we cut off the DLVO po- 
tentials and model the minimum by a parabola 
as described in Refs. [23 HE]- To ensure ex- 
cluded volume of the particles we use a repul- 
sive (Hertzian) potential. Below the resolution 
of the SRD algorithm short range hydrodynam- 
ics is corrected by a lubrication force within the 
MD framework as explained in Refs. [231 I2T1 



For the simulation of a fluid solvent, many 
different simulation methods have been pro- 
posed: direct Navier Stokes solvers [3TJ [321 [33J 
134"] . Stokesian Dynamics (SD) [3J [351 [35] , Ac- 
celerated Stokesian Dynamics (ASD) [37] 135] . 
pair drag simulations [4], Brownian Dynam- 
ics (BD) [T71 EH], Lattice Boltzmann (LB) [fl 
[21 [301 SI], and Stochastic Rotation Dynam- 
ics (SRD) [23 [31 S3J- These mesoscopic 
fluid simulation methods have in common that 
they make certain approximations to reduce the 
computational effort. Some of them include 
thermal noise intrinsically, or it can be included 
consistently. They scale differently with the 
number of embedded particles and the complex- 
ity of the algorithm differs largely. In partic- 
ular, there are big differences in the concepts 
how to couple the suspended particles to the 
surrounding fluid. 

We apply the Stochastic Rotation Dynamics 
method (SRD) introduced by Malevanets and 
Kapral [331 EH]. It intrinsically contains fluctua- 
tions, is easy to implement, and has been shown 
to be well suitable for simulations of colloidal 
and polymer suspensions [24 ] l27 l 142 ] |43 [ 146 ] 147 ] 
HH] , The method is also known as "Real-coded 
Lattice Gas" [32] or as "Multi-Particle-Collision 
Dynamics" (MPCD) [13 [50]. It is based on so- 
called fluid particles with continuous positions 
and velocities. A streaming step and an inter- 
action step are performed alternately. In the 
streaming step, each particle i is moved accord- 
ing to 



r ! (t + r) = r l (t) + rv ! (i), 
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where Yi(t) denotes the position of the parti- 
cle i at time t and r is the time step. In the 
interaction step the fluid particles are sorted 
into cubic cells of a regular lattice and only 
the particles within the same cell interact with 
each other according to an artificial interac- 
tion. The interaction step is designed to ex- 
change momentum among the particles, but at 
the same time to conserve total energy and to- 
tal momentum within each cell, and to be very 
simple, i.e., computationally cheap. Each cell 
j is treated independently: first, the mean ve- 
locity Uj(t') = J2i=i Vi(t) in cell j is 
calculated. Nj(t') is the number of fluid parti- 
cles contained in cell j at time t' = t + r. Then, 
the velocities of each fluid particle in this cell 
are rotated according to 

v i {t + T) = u j (t')+n j (t')- [Vi (t) - Uj (<')]■ (4) 

f2j(<') is a rotation matrix, which is indepen- 
dently chosen at random for each time step and 
each cell. We use rotations about one of the co- 
ordinate axes by an angle ±a, with a fixed. 
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The coordinate axis as well as the sign of the 
rotation are chosen at random, resulting in 6 
possible rotation matrices. However, there is a 
great freedom to choose the rotation matrices. 
Any set of rotation matrices satisfying the de- 
tailed balance for the space of velocity vectors 
could be used here. To remove anomalies intro- 
duced by the regular grid, one can either choose 
the mean free path to be sufficiently large or 
shift the whole grid by a random vector once 
per SRD time step [ST1I52]. 

To couple the SRD and the MD simulation, 
basically three different methods have been in- 
troduced in the literature. Inoue et al. pro- 
posed a way to implement no slip boundary 
conditions on the particle surface [42,. To 
achieve full slip boundary conditions, Lennard- 
Jones potentials can be applied for the inter- 
action between the fluid particles and the col- 
loidal particles (45] [53] . A more coarse grained 
method was originally designed to couple the 
monomers of a polymer chain to the fluid [MJ 
[55] . but in our previous work [22J, [24] [27] we 
have demonstrated that it can also be applied 
to colloids, as long as no detailed spacial res- 
olution of the hydrodynamics is required. We 
use this coupling method in our simulations and 
describe it shortly in the following. 

For the coupling of the colloidal particles to 
the fluid, they are sorted into the SRD cells 
and included in the SRD interaction step. The 
stochastic rotation is performed in the momen- 
tum space instead of the velocity space to take 
into account the difference of inertia between 
light fluid and heavy colloid particles. We have 
described the simulation method in more detail 
in Refs. OEllEH]. 

In the present work we report on several stud- 
ies about pressure filtration, cluster formation 
in a steered simulation, and sedimentation. All 
these studies are proofs of principle, and there- 
fore small simulations are performed. Typical 
parameters for these simulations are: 60 fluid 
particles per box, with 0.296 fim extension in 
each direction. The volume of the box is set 
to be the same as the volume of a colloidal 
particle with diameter 0.37 /im. The system 
size is usually 30 to 60 boxes in each dimen- 
sion. By matching the diffusion constant, den- 
sity and viscosity to a real suspension (23], the 
scaling scheme we have presented in Ref. [27] 
yields a scaled temperature in the simulation of 
14 mK and a scaled viscosity of 4.78 • 10~ 8 Pas. 
To preserve the correct dynamics, characterized 
by the dimcnsionless numbers (Re, Pe, Kn,. . . ) 
one has to rescale the potentials and all driving 
forces in the MD scheme by the same scaling 
factor, as well. 



Let us now shortly sketch the technical re- 
alization of the steering interface in our sim- 
ulation code. The program is an object ori- 
ented code written in C++. Each object 
contains virtual routines saveO and loadO, 
which write the whole object data to or load 
it from a buffer, respectively. The buffer con- 
tains a plain text description of all variables 
contained in the object with their values, sim- 
ilar to the C++ source code one would write 
to initialize an instance of the class. The whole 
work flow of a simulation, including the actual 
MD loop as well as data input/output tasks 
are described in this manner using specialized 
workst ep-classes |59] . each providing a special- 
ized work () -routine which performs different 
actions depending on the actual class type of 
the respective object. One of the workstep- 
classes is designed to change a specified ob- 
ject by using its saveO and loadO routine. 
First, the current values are stored to a tem- 
porary buffer, then one or several variables 
may be overwritten by new values, and finally 
the object to be modified is loaded again from 
the temporary buffer. The description of the 
changes may be read from standard input or 
from a file, similar to the simulation setup, 
which is read during the initialization of any 
simulation, workst ep-class objects can also be 
included into the work flow at later times, or 
may even be disconnected from the usual work 
flow and bound to standard UNIX system sig- 
nals. By default a workstep writing particle 
positions to standard output and a workstep to 
change objects getting its buffer from standard 
input are bound to the SIGTTOU ("terminal out- 
put") and SIGTTIN ("terminal input") signal, 
respectively. This allows to embed the simula- 
tion program into a framework of shell scripts 
which generate the appropriate input, redirect 
the output to a visualization tool and send the 
signals according to the user interaction. This 
can be realized in a client-server fashion, even 
on different hosts and platforms using appro- 
priate scripts and TCP/IP connections. 



III. RESULTS 

After the more technical paragraph in the 
previous section we now turn to a discussion of 
the advantages of such a steering approach for 
computer simulations and highlight some pit- 
falls resulting from the physical background in 
the context of steering. This section is subdi- 
vided into several subsections each of them fo- 
cusing on a particular example to illustrate the 
steering approach in practice. 
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Figure 1: Local density expressed in terms of vol- 
ume fraction of the solid particles constituting the 
filter cake depending on the pH-value. The ionic 
strength is kept constant at / = 3mmol/l, and the 
volume fraction at $ = 5%. The plots for different 
pH-values are shifted against each other vertically 
by 0.2 for better visibility. The shape of the profile 
differs, since the spacing of the particles depends 
on their interactions. 



A. Pressure filtration 

In this case a filtration of a suspension is sim- 
ulated. The suspended particles cannot pass 
the filter, whereas the fluid passes through the 
filter without resistance in an idealized filtra- 
tion process. The suspended particles agglom- 
erate in front of the filter and form a so-called 
filter cake. Since the dynamics of the parti- 
cles is not only governed by the hydrodynamics 
of the fluid, but also by their DLVO interac- 
tion, the density and the structure of the filter 
cake will depend on the pH-value and the ionic 
strength. 

In our simulation we drive the fluid by an 
applied force pointing downwards acting in a 
small region close to the upper boundary of 
the system. For the fluid we apply fully peri- 
odic boundary conditions, whereas the bound- 
aries for the suspended particles are closed in 
z-direction. In this way, the fluid is forced to 
stream in vertical direction and to drag the par- 
ticles to the bottom of the system. In test sim- 
ulations we have chosen constant ionic strength 
/ = 3mmol/l and a constant overall volume 
fraction of <f> = 5%, and simulated the filtra- 
tion process for several pH-values. 

Depending on the interactions, different in- 
ternal structures of the filter cake are formed 
when the density of the particles increases at 
the bottom of the system. To gain a good 
statistics for the structure large simulations are 
needed. However, in a filtration process, usu- 



ally low initial densities are used, so that very 
large simulation volumes arc needed to obtain 
large particle numbers for good statistics. On 
the other hand, large simulation volumes also 
imply large data files and long transient times 
until the filter cake is formed and a considerable 
part of the particles has reached the bottom of 
the system. 

This is a typical problem, in which one would 
like to observe the running simulation occasion- 
ally to check if the filter cake is already formed 
and when the structure does not change any- 
more. One would like to measure pressure pro- 
files, local streaming velocities, or simply the 
final density profile of the filter cake. Here, 
steering would mean to initiate data acquisi- 
tion by user interaction. This is of advantage if 
an action, like starting to average certain data, 
should be initiated when conditions arc fulfilled, 
which are difficult to check automatically. Since 
the density and structure of the filter cake is a 
priori unknown, an automated check is difficult 
to implement. 

Additionally, the upper boundary of the fil- 
ter cake is diffuse and depends on the condi- 
tions (pH-value and ionic strength) as well. In 
Fig. [I] the dependence of the density profiles of 
filter cakes obtained in simulations for different 
pH- values are shown. With increasing pH- value 
the height of the filter cake increases and the 
density profile changes slightly. Some voids in 
the filter cake diminish the local density if com- 
pared to a dense packing. At the bottom of 
the filter cake, layers of particles are present, 
whereas in the upper regions the structure be- 
comes more irregular. However, the structure 
of the sediment and its height will also depend 
on the pressure exerted on the fluid and proba- 
bly on the initial volume fraction. Larger sim- 
ulations than this proof of principle are needed 
to be able to quantify the pH- value-dependence 
of the height of the filter cake. However, sim- 
ulations without driving force on the fluid, can 
also help to understand this system. Therefore, 
we explore the parameter space in the absence 
of a driving force in the following subsection. 



B. Observation of cluster formation 

Since several parameters, namely the vol- 
ume fraction, pH-valuc, ionic strength, exter- 
nal driving forces influence the system, a first 
step to understand its behavior is to exam- 
ine it in the absence of external forces. Still, 
two parameters-the pH-value and the ionic 
strength-govern the dynamics of the system. If 
the potentials are attractive, cluster formation 
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Figure 2: Memory effects in a steered simulation: the grayscale denotes the coordination number. Particles 
closer than 2.4 radii are considered to be neighbors. Dark particles in this sense are highly coordinated, 
bright ones have low coordination numbers, a) right after initialization, homogeneous distribution of the 
particles, b) cluster formation after tuning the potentials to be attractive, c) still inhomogeneous regions 
after steering the simulation out of the clustered region, d) the system is homogeneous again when diffusion 
has counterbalanced the density inhomogeneities. The steering path in the parameter space as well as the 
points a)-d) thereon are illustrated in Fig.[3] 



can be observed [57], provided the attraction 
is strong enough compared to thermal fluctu- 
ations. In our previous work |221 128] we have 
explored this part of the parameter space for 
some distinct volume fractions by running nu- 
merous individual simulations. However, with 
the newly implemented steering approach one 
can detect boundaries between different regions 
of the stability diagram faster and more sensi- 
tively. 

In Fig. [2] we illustrate the effects of chang- 
ing the properties of the interaction potentials. 
All the snapshots belong to the same simula- 
tion and represent different points of the tra- 
jectory the system is steered through param- 
eter space: a) The initial situation with repul- 
sive potentials, shortly after initialization of the 
system, b) after crossing the boundary to the 
clustered region, c) shortly after steering back 
to the suspended regime, where still some in- 
homogeneities in the system can be observed, 
and finally d) after some time when diffusion 



has restored a homogeneous system. This time 
depends on the distances between the clusters 
and on the characteristic diffusion time of the 
particles. The steering path in the configura- 
tion space is shown in Fig.[3l The points at 
which the snapshots of FigT|2] are taken, are 
marked by black dots. In the insets in Fig. [3] 
the total potential is plotted for several cases. 
During the simulation we change the pH-value 
and the ionic strength and we drive the system 
along a closed trajectory in parameter space. 
We start at low ionic strength / = lmmol/1 
and pH — 4 (potential shown) and first increase 
the pH- value (potential shown for 1=1 mmol/1 
and pH = 7), then increase the ionic strength, 
then reduce the pH-value again while further 
increasing the ionic strength (potential shown 
for / = 15 mmol/1 and pH = 4) and finally de- 
crease the ionic strength to return to the start- 
ing point. We choose the trajectory such that 
the barrier between the primary and the more 
shallow, secondary minimum of the DLVO po- 
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tential does not go below 5 k^T. Since the sec- 
ondary minimum becomes deeper during the 
simulation cluster formation can be observed. 
However, the barrier between the primary and 
the secondary minimum of the DLVO poten- 
tials prevents irreversible aggregation in the 
primary minimum. A case in which the bar- 
rier vanishes and irreversible aggregation would 
appear is also shown in the upper right inset 
(I = 15mmol/l and pH — 7). In any case, 
in our simulation all clusters dissolve again af- 
ter steering back to a lower pH- value and lower 
ionic strength (see Fig.[2]d)). This confirms 
that we have reversible clustering in the sec- 
ondary minimum of the DLVO potentials. 

One of the advantages of a steered simula- 
tion is, that one can more accurately detect the 
onset of the cluster formation than by start- 
ing individual simulations for different condi- 
tions. One might predict the boundary by eval- 
uating the depth of the secondary minimum in 
the potential, but since also the hydrodynamic 
damping forces and an eventually applied driv- 
ing pressure influence this process the predic- 
tion of the boundary can become more compli- 
cated. We have experienced similar difficulties 
in the context of shear flow simulations of the 
same colloidal suspension [22l 124] . 

A quantity to support the visualization of 
clusters is the coordination number. We con- 
sider two particles whose centers are closer than 
2.4 particle radii as neighbors and define the 
coordination number of each particle as the 
number of its neighbors in the sense just de- 
scribed. The particles in Fig. [2] are drawn in 
a grayscale corresponding to their coordination 
number, where dark particles are highly co- 
ordinated particles and bright ones are those 
with a low coordination number. At the begin- 
ning of the simulation (top left) the particles 
are distributed homogeneously over the whole 
simulation volume. After steering the simula- 
tion in the parameter space into the region of 
attractive potentials, cluster formation can be 
observed (top right). When the potentials are 
made less attractive again, the clusters dissolve 
(bottom left). However, the particles prefer- 
ably stay at their positions and thus, the in- 
homogeneities in the system do not disappear 
immediately, which can be seen also in the dark 
color denoting high coordination numbers. Af- 
ter some time, diffusion restores the homogene- 
ity of the system (bottom right). However, if 
the pH-value or the ionic strength is increased 
too much, so that the barrier between the pri- 
mary and the secondary minimum of the DLVO 
potentials vanishes, the clustering process be- 
comes irreversible, which we avoided by the 



choice of the steering path (Fig. [3}. 

In this example the simulation remembers 
the interaction by the user. Even if the steering 
path is selected carefully inside the reversible 
range, the simulated suspension needs a char- 
acteristic time to relax after changing the in- 
teractions. The interaction by the user can be 
seen as a perturbation in the physical sense and 
the system needs time to adopt to the new sit- 
uation. In these cases special care is advised 
when steering a simulation. We illustrate this 
by another example in the following subsection. 



C. sedimentation: hydrodynamic 
interaction 

Not only the particle positions and velocities, 
but also hydrodynamic interactions can influ- 
ence the behavior of the system. Therefore, in- 
homogeneous particle distributions induced by 
steering the simulation first into a clustered re- 
gion of the parameter space before simulating 
a different situation, might influence the result: 
we show in the following paragraph that a sed- 
imentation process can happen faster, if clus- 
ter formation by attractive interactions creates 
density inhomogencitics in the system. 

The friction force a single particle feels in a 
sedimentation process can be calculated analyt- 
ically, but it differs if other particles are present. 
In Ref. [57] we have confirmed that the sedimen- 
tation velocity in our simulation depends on the 
volume fraction, as it is well-known from sed- 
imentation theory |58j . But even at constant 
volume fraction the sedimentation velocity can 
be different in different configurations. If the 
particles form clusters, the whole cluster settles 
down and the fluid streams around the whole 
cluster. The resistance is much less compared 
to the case when the particles are distributed 
homogeneously and the fluid streams around 
each of them separately (compare Ref. [55] and 
references therein). In Fig.|4] we have plotted 
the sedimentation velocity evaluated in several 
simulations, which only differ by the potentials. 
The volume fraction is the same for all of them. 
However, as one can see in the figure, the sedi- 
mentation velocity is smaller for low pH- values. 

This can be explained as follows: Since for in- 
creasing pH- value the potentials become attrac- 
tive, the particles form clusters and settle down 
faster. This effect is purely due to the hydro- 
dynamics of the system. Since the streaming 
field depends on all particle positions and veloc- 
ities not only at a given time, but also on their 
history, steering in this context may be very 
dangerous. When long range interactions, like 
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Figure 3: Steering path of the simulation: we start at the bottom left corner at low pH-value and low 
ionic strength and steer the simulation along the dashed path through parameter space. The region of 
reversible cluster formation is reached, but the barrier between primary and secondary minimum of the 
DLVO potentials prevents irreversible clustering: when reaching the starting point again, all clusters have 
vanished. The boundary of the clustering region is sketched: Reversible clustering is observed when the 
depth of secondary DLVO minimum is in the order of ksT and irreversible clustering appears when the 
height of the barrier between the two DLVO minima is of the order of k^T or less. The labeled black dots 
mark the positions at which the snapshots of Fig. [2] are taken. The insets show the total potential at some 
selected extreme cases (see text). 
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Figure 4: Sedimentation velocity of a $ = 5% sus- 
pension in a closed vessel. The sedimentation veloc- 
ity of the particles in the upper part of the system 
is averaged over several time steps. The particles 
at the bottom of the system are not encountered. 
The velocity depends on the particle-particle inter- 
action. Clusters settle down faster than individual 
particles. Symbols denote simulation results, the 
line is a guide to the eye. 

hydrodynamics, which also depends on the his- 
tory of the system, is important, steering should 
be avoided. In this case, to obtain reasonable 
results, one has to start again from an initial 



configuration when changing the interactions. 



IV. CONCLUSION 

Bearing in mind the underlying physics, one 
can roughly estimate the limits of steered sim- 
ulations. Let us summarize the points we have 
brought forward in this paper. First of all, 
steering can save computing time when aiming 
at a rough understanding of the influences on a 
system, especially when looking for the "inter- 
esting points" in parameter space, e.g., when 
searching for transition lines between different 
phases. We have shown this by the onset of 
cluster formation due to modifying the interac- 
tion potentials in a model suspension. Steering 
can be used as well to start data acquisition af- 
ter a transient at the beginning of a simulation 
or to adjust the frequency for the data acquisi- 
tion according to the current state of the sim- 
ulation. This is especially interesting for large 
simulations with long transient times, as illus- 
trated by the example of filter flow. We also ex- 
plained why in some cases steering may be seen 
as a perturbation of the system, especially when 
the interaction potentials are changed. Special 
care is advised when dealing with non ergodic 
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systems, memory effects, or long range inter- 
actions, like the role of hydrodynamics in sedi- 
mentation in our last example. 
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